Phase diagram of a two- component Fermi gas with resonant interactions 
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The pairing of fermions is at the heart of superconductivity and superfluidity. The recent ex- 
perimental realization of strongly interacting atomic Fermi gases has opened a new, controllable 
way to study novel forms of pairing and superfluidity. A major controversial issue has been the 
stability of superfluidity against an imbalance between the two spin components when the fermions 
interact resonantly. Here we present the phase diagram of a spin-polarized Fermi gas of ^Li atoms 
at unitarity, mapping out the superfluid phase versus temperature and density imbalance. Using 
tomographic techniques, we reveal spatial discontinuities in the spin polarization, the signature of 
a first-order superfluid-to-normal phase transition, which disappears at a tricritical point where the 
nature of the phase transition changes from first-order to second-order. At zero temperature, there 
is a quantum phase transition from a fully-paired superfluid to a partially-polarized normal gas. 
These observations and the implementation of an in situ ideal gas thermometer provide quantitative 
tests of theoretical calculations on the stability of resonant superfluidity. 



Superfluidity and superconductivity of fermions are 
based on the formation of fermion pairs. The stabil- 
ity of these pairs determines the robustness of the su- 
perfluid state, and the quest for superconductors with 
high critical temperature is a search for systems with 
strong pairing mechanisms. Ultracold atomic Fermi gases 
present a highly controllable model system for study- 
ing strongly interacting fermions [1 . Tunable interac- 
tions utilizing Feshbach collisional resonances and con- 
trol of population or mass imbalance among the spin 
components provide unique opportunities to investigate 
the stability of pairing [2} El S], and possibly to search 
for new exotic forms of superfluidity [5] [6]. The case 
of unitarity, when the two spin components resonantly 
interact and the behavior of the system becomes inde- 
pendent of the nature of the interactions, has become a 
benchmark for experimental and theoretical studies over 
the last few years. However, there is an ongoing debate 
about the stability of resonant superfluidity, reflected in 
major discrepancies in predicted transition temperatures 
for the balanced spin mixture [71 [8| [9], and an even 
more dramatic discrepancy for the critical imbalance 
of the two spin components, called the Chandrasekhar- 
Clogston (CC) limit of superfluidity [2 , 3 . Recent Quan- 
tum Monte-Carlo (QMC) calculations predicted that su- 
perfluidity would be quenched by a density imbalance 
around 40% [10^, whereas other studies predicted a crit- 
ical imbalance^bove 90% [H Ulini El [IH [16]. Our 
earlier work [iTl [181 HI] suggested the lower limit but 
other experiments ^ were interpreted to be consis- 
tent with the absence of the CC limit. This huge dis- 
crepancy reveals that even qualitative aspects, such as 
the role of interactions in the normal phase, are still con- 
troversial. The lack of reliable thermometry for strongly 
interacting systems limits the full interpretations of ex- 
perimental results. 

Here we resolve this long standing debate by present- 



ing the phase diagram of a spin-polarized Fermi gas 
at unitarity. We observe that the normal-to-superfluid 
phase transition changes its nature. At low tempera- 
ture, the phase transition occurs with a jump in the 
spin polarization as the imbalance increases, which we 
interpret as a first-order phase transition. The local 
spin polarization or local density imbalance is defined 
as cr = (n| — ni)/{n^ +^1)7 where t and | refer to the 
two spin components with densities n^^^. At high tem- 
perature, the phase transition is smooth and therefore 
of second-order. The two regimes are connected by a 
tricritical point [41 |23l for which we estimate the po- 
sition as ((Jtc, Ttc/Tp|) ^ (0.2,0.07), where /cb^ft = 
fi^(67r^n|)^/^/2m is the Fermi energy of the majority 
component of density (ks is the Boltzmann constant, 
ft is the Planck constant divided by 27r and m is the 
atomic mass of ^Li). Our low-temperature results con- 
firm a zero-temperature quantum phase transition at a 
critical polarization (Jco ^ 0.36. 

This work required the introduction of several novel 
techniques. A tomographic reconstruction of local Fermi 
temperatures and spin polarization allowed us to ob- 
tain the phase diagram for the homogeneous system, 
no longer affected by the inhomogeneous density of the 
trapped samples. Furthermore, absolute temperatures 
were obtained using in situ thermometry applied to the 
non-interacting fully-polarized Fermi gas in the outer 
part of the trapped samples, an ideal thermometer with 
exactly-known thermal properties. In contrast to previ- 
ous works [18, 22 , this is a direct measurement without 
any approximations. 

Spatial structure of a trapped Fermi mixture 

Our experiments are carried out in a trapping potential 
V{r). The local chemical potential of each spin compo- 
nent is given as /it,|(r) = /i|o,io ~ where /i|o,io 
are the global chemical potentials. When /i^o 7^ Mio due 
to imbalanced populations, the chemical potential ratio 
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FIG. 1: Schematic of spatial structure of a strongly interact- 
ing Fermi gas in a harmonic trap, (a) A two-component (spin 
t and I) Fermi mixture is confined in an external potential 
V(r) oc with the chemical potentials of each spin compo- 
nents /ito,io (^Mi is the shift for the spin | component due to 
interactions), (b) Density distributions of the majority com- 
ponent n|(r) (red) and the minority component n|(r) (blue), 
(c) Spin polarization a{r) = (n^ — ni)/{n^ + ^i)- At zero 
temperature, the trapped Fermi mixture has radially a three- 
layer structure: (I) The core region (0 < r < Rc) of a fully 
paired superfluid with equal densities of the two components, 
(II) the intermediate region (Rc < r < i?|) of a partially po- 
larized normal gas, and (III) the outer region {R^ < r < R^) 
of a fully polarized normal gas. The critical polarization ac 
(as) is defined as the minimum (maximum) spin polarization 
of the normal (superfluid) region. The non-interacting case is 
shown in insets, a.u., arbitrary units. 
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FIG. 2: Double in situ phase-contrast imaging of a trapped 
Fermi mixture. Two phase-contrast images of one sample 
were taken using different probe frequencies of the imaging 
beam, measuring (a) the density difference ridi — —rii and 

(b) the weighted density difference nd2 = 0.76n| — 1.43n|, 
respectively. The phase- contrast images show the 2D dis- 
tribution of the column density difference, ndi,2{x,z) = 
J dy ndi,2{r) where the integral describes the line-of-sight in- 
tegration. The field of view for each image is 150 /imx820 fim. 

(c) The distributions of the column density difference hdi 
(black) and nd2 (red) along the central line (the dashed 
lines in (a) and (b)). The profiles of the integrated lin- 
ear density difference, (d) ndi,z = f dx ndi{x,z) and (e) 
ndi,x = J dz ndi{x,z) show the identical flattop feature ex- 
cept scaling. The aspect ratio of the trapping potential was 
A = 6.15, the majority atom number was = 5.9(5) x 10^, 
the population imbalance was S = 44(4)%, and the relative 
temperature was T' — T /Tfq — 0.03(1) (see text for defini- 
tions). 



77(r) = iJiJiJi^ spatially varies over the trapped sample 
and thus, under the local density approximation (LDA) 
the trapped inhomogeneous sample is represented by a 
line in phase diagrams of the homogeneous system. Fig- 
ure [T] illustrates the spatial structure of a strongly in- 
teracting Fermi mixture in a harmonic trap. In the in- 
ner region, where r] is closer to unity, a superfluid with 
zero (small) spin polarization will form at zero (low) 
temperatures, having a sharp phase boundary against 
the partially-polarized normal gas in the outer region. 
The spin polarization shows a discontinuity at the phase 
boundary r = Rc, a signature of the phase separation of 
a superfluid and a normal gas [24 . The critical polariza- 
tion cTc = lim^^^+ cr(r) {as = lim^^^- cr(r)) represents 
the minimum (maximum) spin polarization for a sta- 
ble normal (superfluid) gas. At higher temperature, the 
discontinuity in the density imbalance disappears. The 
main result of this paper is the observation and quanti- 
tative analysis of such density profiles. Since we have no 
experimental evidence, we are not discussing the exotic 



partially-polarized phases [29^ which could exist only in 
the transition layer between the superfluid core and the 
normal outer region. 

We prepared a variable spin mixture of the two lowest 
hyperfine states of ^Li atoms, labeled as | j) and | |), 
at a magnetic field of 833 G. A broad Feshbach reso- 
nance is located at 834 G and the interactions between 
the two spin states are resonantly enhanced. Our sample 
was confined in a 3D harmonic trap with cylindrical sym- 
metry. The in situ density distributions of the majority 
(spin I) and minority (spin |) components were deter- 
mined using a phase-contrast imaging technique. Since 
the trapped sample was observed to have an elliptical 
shell structure of the same aspect ratio A = fp/fz as 
the trapping potential over our entire temperature range, 
where fp {fz) is the oscillation frequency in the transverse 
(axial) direction (Fig. [2|, we obtained the low- noise pro- 
files n(r) by averaging the column density distribution 
along the equipotenital line (defined as X^x^ + = 
for a given radial position r). The region for averaging 
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FIG. 3: Density profiles of trapped Fermi mixtures with imbalanced populations. The first row (a-d) shows the averaged 
column density profiles for various temperatures (red: majority, blue: minority, black: difference). The majority radius 
was determined from the outer region (r > i?|: the radius of the minority cloud) of the majority profiles using a fit to a 
zero-temperature Thomas- Fermi (TF) distribution (black dotted lines). The column densities are normalized by the central 
value of the fitted TF distribution. The second row (e-h) and the third row (i-1) show the reconstructed 3D profiles and the 
spin polarization profiles a{r) corresponding to the profiles in a-d. The core radius Rc was determined as the peak (and/or 
kink) position in the column density difference (only for a-c), indicated by the up arrows and the dashed lines. The two spin 
polarizations ctc at r = Rc and cts at r = i?c — 0.05i?| are marked by the right and left arrows, respectively. T\ ac, Rc/R^, R^ 
(in /xm), N^, 6 (in %) and A were respectively: (a, e, i) 0.03(1), 0.34, 0.43, 385, 5.9(5) x 10^ 44(4), 6.15; (b, f, j) 0.05(2), 0.24, 
0.39, 416, 1.0(1) X 10^ 48(4), 6.5; (c, g, k) 0.07(1), 0.21, 0.29, 443, 1.2(2) x 10^ 54(4), 6.5; (d, h, 1) 0.10(1), not determined, 
not determined {ar=o = 0.15 and condensate fraction = 2(1)%), 398, 5.3(4) x 10^, 54(4), 7.7. 



was restricted depending on the type of analysis. The 
density profiles n(r) were determined from the 3D recon- 
struction using the inverse Abel transformation of the 
column densities n(r) [25]. Deviations from the trap as- 
pect ratio were only found for the outer thermal wings 
and will be discussed below. More detailed description 
of the experimental procedure and the image process- 
ing is provided in appendix. Most of our measurements 
were performed at a total population imbalance S ^ 50%, 
where S = {N^ - 7V|)/(7V| + N^) refers to the total atom 
numbers in the sample, and of the spin | and | 
components, respectively. 

Figure [3] displays the radial profiles of the densities 
n|(r) and n^(r) of each component and the correspond- 
ing spin polarization cr(r) for various temperatures. The 
discontinuity in the spin polarization, clearly shown at 
very low temperature, demonstrates the phase separa- 
tion of the inner superfluid of low polarization and the 



outer normal gas of high polarization. At low tem- 
perature the core radius Rc is determined as the kink 
(and/or peak) position in the column density difference 
profile. At high temperature (but still in the superfluid 
regime), the discontinuity in a{r) disappears. At our 
lowest temperature, the radii of the minority cloud and 
the core region were measured as R^ = 0.73{1)R^ and 
Rc = 0.430(3)i?| (at S = 44(4)%), respectively, and agree 
with recent theoretical calculations [TOl |29j within the ex- 
perimental uncertainties due to the determination of 6. 
Here, R^ is the radius of the majority cloud, measured 
as R^ ^ 0.95Rtf {Rtf is the Thomas-Fermi (TF) ra- 
dius of a non-interacting Fermi gas with the same atom 
number of A^|), reflecting the attraction between the two 
components. 

Thermometry of a strongly interacting Fermi gas 

In most cold atom experiments, temperature is deter- 
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FIG. 4: Temperature determination using in situ density pro- 
files. The relative temperature T' = T/Tfq (see text for defi- 
nition) was determined from the outer region {r > Ri) of the 
averaged column density difference profile (black line) fitted 
to a finite temperature Fermi-Dirac distribution (red line). 
The radius of the minority cloud Ri was determined from a 
fit of the wing profile of the minority component (black dashed 
line) to a zero temperature TF distribution (red dashed line) . 
(a) T' = 0.03(1) and 5 = 44(4)%, and (b) T' = 0.08(1) and 
3 = 46(4)%. 



mined from the density distribution after ballistic ex- 
pansion, reflecting the momentum distribution of ther- 
mal atoms in the cloud. In the case of a population- 
imbalanced Fermi mixture, the outer part of the major- 
ity component, having no spatial overlap with the minor- 
ity component, is a non-interacting Fermi gas in thermal 
equilibrium with the inner part and fulfills the defini- 
tion of an ideal thermometer, namely a substance with 
exactly-understood properties in contact with the sam- 
ple to be characterized. This concept was first introduced 
in ref. [18 , where temperature was determined from the 
majority wing profile after expansion. However, we found 
that during expansion some of the outer majority atoms 
experience collisions with the minority atoms in the in- 
ner region, causing a modification of their kinetic energy 
(see appendix). We avoided this problem by analyzing 
the in situ profiles. The outer part of the averaged col- 
umn density difference profile (r > R^) was fit to a fi- 
nite temperature Fermi-Dirac distribution in a harmonic 
trap (Fig [4| and the relative temperature T' = T/Tpo 
was determined, where UbTpq = fi^(67r^no)^'^^/2m is 
the Fermi energy of the non-interacting Fermi gas which 
has the same density distribution in the outer region as 
the majority cloud (no is the central density of the non- 
interacting Fermi gas at zero temperature). We verified 
that anharmonicity of the trapping potential does not 
affect the fitted temperature (see appendix). 

Phase diagram for a homogeneous system 

The critical lines of the phase diagram of a homoge- 
neous spin-polarized Fermi gas were obtained by de- 



termining the local temperature and spin polarization 
at the phase boundary. The local relative temperature 
T/^^ = T/Tp] was derived from the local density n^{Rc) 
as V{Rc) = T/Tpo x {no/n^{Rc))^^^ . Since we observe 
no jump in the majority density within our resolution, 
Ti?| is well-defined at the phase boundary. The critical 
spin polarizations ac and as were measured as ac = cr{Rc) 
and as = a{Rc — O.ObR^) [37 . The discontinuity in the 
spin polarization profile implies that there is a thermo- 
dynamically unstable window, as < a < ac^ leading to a 
first-order superfluid-to-normal phase transition. As the 
temperature increases, the unstable region reduces with 
(Jc decreasing and as increasing. For high temperature 
when the bimodal feature in the spin polarization pro- 
file disappears, we recorded the condensate fraction as 
an indicator of superfluidity, using the rapid field-ramp 
technique pTl. As the temperature decreases, the con- 
densate fraction gradually increases with a finite central 
polarization [19]. Such a smooth variation of the density 
profile and condensate fraction across the phase transi- 
tion are characteristic of a second-order phase transition. 

The phase diagram is characterized by the three dis- 
tinguished points: the critical temperature Tco for a bal- 
anced mixture, the critical spin polarization aco of a nor- 
mal gas at zero temperature, and the tricritical point 
((Jtc, Ttc) where the nature of the phase transition changes 
from second-order to first-order. Due to the lack of a 
predicted functional form for the phase transition line in 
the a-T plane, we apply a linear fit to the measured crit- 
ical points, suggesting Tco/Tf^ ^ 0.15, aco ~ 0.36 and 
{crtc,Ttc/TF^) ^ (0.20,0.07). The value for aco agrees 
well with the prediction of QMC calculation of 0.39 [10^. 
The extrapolation of the phase diagram to a = is only 
tentative, since the in situ thermometry could not be 
applied to small population imbalances due to the nar- 
rowness of the non-interacting outer region. It is possible 
that the first- and second-order transition lines meet at 
the tricritical point with different slopes, similar to the 
case of liquid ^He-'^He mixtures [26 . 

The zero-temperature phase diagram and in particu- 
lar the value for the CC limit are significantly affected by 
strong interactions in the normal phase. The CC limit 
reflects the energetic competition between a superfluid 
state and a partially-polarized normal state. When the 
chemical potential difference Sja = ja^ — /i| is larger than 
a critical difference 2hc, the normal state is energetically 
favorable and the superfluid state breaks down. In BCS 
theory, valid for weak interactions, the critical difference 
is he = A/a/2 (A is the pairing gap). At unitar- 
ity, QMC studies predict he 1.00(5)A(?:^ 1.2/i) HI] 
with the assumption of a non-interacting normal gas. 
Here, /i = (/i^ + The condition /i^c = /J^ — he < 

implies that = in a non-interacting normal 
gas, and consequently <Jco = 100%, i.e. the absence 
of a partially-polarized normal phase. Mean- field ap- 
proaches ^121 [131 HH [13 [IS! 5 which cannot treat the in- 
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FIG. 5: a-T phase diagram for a homogeneous spin-polarized 
Fermi gas with resonant interactions. The critical polariza- 
tions (Tc (black solid circles and square) and as (gray solid cir- 
cles) are displayed along the local T/Tp^ at the phase bound- 
ary. The yellow area {as < a < ac) represents a thermody- 
namically unstable region, leading to the phase separation. 
Above the tricritical point, the phase transition in the center 
of the cloud was observed by the onset of pair condensation. 
For this, a cloud was evaporatively cooled, until it crossed 
the phase transition on a trajectory almost perpendicular to 
the phase transition line (see appendix). The critical spin 
polarization and temperature were obtained by interpolating 
between points without and with small condensates (black 
solid square). The linear fit to the acS is shown as a guide 
to the eye for the normal-to-superfiuid phase transition line. 
Each data point consists of five independent measurements 
and error bars indicate standard deviation. The blue open 
symbols show theoretical predictions for the critical tempera- 
ture of a homogeneous equal mixture (v- Bulgac et al. [7], A: 
Burovski et al. [8 , O: Haussmann et al. 9 ) and the critical 
polarization at zero temperature (□: Lobo et al. [IQ]). The 
blue solid square is the measured critical temperature from 
Luo et al. |22|, multiplied by with ^ — 0.42 [TI] to obtain 
local T/Tf at the center. Finite temperature correction may 
increase the effective value of ^. 



teractions in the normal phase in an accurate way, also 
predict a high critical imbalance (Jco > 90%. Strong in- 
teractions between the atoms in the normal phase, how- 
ever, have been observed through the compressed shape 
of the minority cloud [18] and the shift in the RF excita- 
tion spectrum [27j. The data in Figure [s] clearly establish 
a zero-temperature CC limit for cjco in the range of 30% 
to 40%. 

The density profiles at our lowest temperature pro- 
vide quantitative information on the zero-temperature 
thermodynamics [28l [29]. At zero temperature, the 
global chemical potential of a fully-paired superfluid in 
the core is given as jiso = £>^f = ^^^(67r^nso)^^^/2m 



where £i? is the local Fermi energy and Ugo is the ma- 
jority (or minority) density at the center, whereas /i^o = 
/i^(67r^no)^/^/2m and /i^o = ^oMto- From the thermo- 
dynamic equilibrium condition /i^o = (mto + Mio)/2, we 
obtain the chemical potential ratio as 

In our coldest sample {5 ^ 45%), the normalized cen- 
tral density and the radii for the phase boundary and 
the minority cloud were measured to be Uso/no = 
1.72(4), i?c/i^T = 0.430(3), and Ri/R^ = 0.728(8), re- 
spectively, yielding r]c = vi^c) ^ 0-03 and r]^ = r]{Ri) ^ 
-0.69 with ^ = 0.42 [11 . Furthermore, the critical dif- 
ference is given as he/ /i = (1 — r]c) /(I + Vc) = 0.95. Since 
theory clearly predicts /i < A O |TT], we have he < A. 
If he were larger than A, polarized quasi-particles would 
have negative energies and form already at zero temper- 
ature. Therefore, up to our observed value of he, the 
fully-paired superfluid state is stable, and a polarized su- 
perfluid exists only at finite temperature. 

The interface between two immiscible fluids involves a 
surface energy, leading to at least a small violation of the 
LDA. However, the observed sharp interface along the 
an equipotential line and the flattop structure of the lin- 
ear density difference profiles (Fig. [2]i and e) imply that 
corrections to the LDA are smaller than the resolution 
of our experiment. These observations are inconsistent 
with the interpretations given for the experimental re- 
sults reported in ref. ^20l [21] , where it has been shown 
that highly-elongated small samples are deformed by sur- 
face tension [30", "SP . The scaling of those surface effects 
to our parameters predicted a deviation of the aspect ra- 
tio of the superfluid core of ^ 15% from the trap aspect 
ratio [31 , whereas we observe this deviation to be smaller 
than 2%. Note that surface tension would add energy in 
the phase-separated superfluid regime and would shift 
the CC limit to smaller values. Ref. [20[ [21] concluded 
that the CC limit should be Seo > 95% which is ruled 
out by our observations. We are not aware of any sug- 
gested effect which can reconcile the data of ref. [20l [21] 
with our phase diagram for a resonant superfluid. To in- 
dentify this finite size effect and to fully understand the 
nature of the normal state [27 are still open questions 
for imbalanced Fermi gases. 

Conclusions 

We have established the phase diagram of a homogeneous 
spin-polarized Fermi gas with resonant interactions in the 
(j-T plane. This includes the identification of a tricritical 
point where the critical lines for first-order and second- 
order phase transitions meet, and the final confirmation 
of a zero-temperature quantum phase transition, the CC 
limit of superfluidity, for a gas at unitarity. So far, pre- 
dicted exotic superfluid states such as the breached-pair 
state in a stronger coupling regime ( "BEC side" ) [I3l [32] 
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and the FFLO state in a weaker coupling regime ("BCS 
side") [H m [331 El EH |36] have not been observed, 
but the novel methods of tomography and thermometry 
will be important tools in the search for those states. 

We thank M. W. Zwierlein and A. Keshet for a critical 
reading of the manuscript. This work was supported by 
NSF and ONR. 



APPENDIX 
Experimental details 

The experimental procedure has been described in pre- 
vious publications [ITl [181 HS] • A degenerate Fermi gas 
of ^Li atoms was first prepared in an optical trap, using 
laser cooling and sympathetic cooling with ^^Na atoms. 
A variable spin mixture of the two lowest hyperfine states 
I t) and I I) (corresponding to the \F = 1/2, = 1/2) 
and \F = 1/2, m^? = —1/2) states at low magnetic field) 
was created at a magnetic field B = 885 G. The final 
evaporative cooling by lowering the trap depth and all 
measurements were performed at B = 833 G. The tem- 
perature of the cloud was controlled by the lowest value 
of the trap depth in the evaporative cooling process. The 
axial trap frequency was fz = 23 Hz. The two transverse 
trap frequencies are equal within less than 2%. 

The optical signal in the phase-contrast imaging is pro- 
portional to the net phase shift of the imaging beam pass- 
ing through a Fermi mixture, i.e. c^n^ — ciui (x / {u — 
v^) -\-ni/{u — i^|), where u is the probe frequency of the 
imaging beam, and z/^ and z/^ are the resonance frequen- 
cies of the optical transition for the states | t) and | |), re- 
spectively. When the probe beam is tuned to the middle 
of the two transitions, i.e. u = uq = {u^ the opti- 

cal signal reflects the density difference Ud = n^—ni with 
C| = C|. In our experiment, two phase-contrast images of 
the same sample were taken consecutively with different 
probe frequencies, ui and 1^2- The two images record the 
density difference Udi = — and the weighted den- 
sity difference nd2 = a^n^ — <^|^|. J^i was determined 
by zeroing the optical signal with an equal mixture and 
a^^l was determined by the signal ratio between the first 
and the second image for a highly imbalanced Fermi mix- 
ture with 1^1 > 95% (an almost fully polarized gas). Fi- 
nally, we obtained = {aiUdi — nd2)/{(^i — a^) and 



ni = {a^Udi — ^d2)/(<^i — a^). The difference between 
ui and 1^2 was chosen to lie between 8 and 13 MHz. The 
time interval between the two images was 10 /is, and the 
pulse duration of each probe beam was 15 /is. Because 
the probe beam was off-resonant, no heating effect of the 
first pulse was observed in the second image. 

Low-noise profiles were obtained by averaging the col- 
umn density distribution of phase-contrast images along 
elliptical equipotential lines (A^x^ + 2:^ = r^). For the 
measurement of the critical spin polarization, the av- 
eraging region was restricted to \x\ < 12 /im in or- 
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FIG. 6: Trajectory of the center of a cloud in the phase dia- 
gram during the cooling process. Above the tricritical point, 
the normal-to-superfluid phase transition was observed by the 
onset of pair condensation in the evaporative cooling process. 
The local spin polarization and temperature at the center of 
the cloud was measured (black solid (open) circles with (with- 
out) condensate fraction) and the critical point was obtained 
by linearly interpolating with the condensate fraction. The 
dashed-dot line shows a guide line for the trajectory of the 
cloud center. The population imbalance of the sample was 
6 ~ 55%. A non-interacting mixture with this imbalance has 
a spin polarization a ^ 30% at the center at zero temperature. 



der to preserve the sharp features at the phase bound- 
ary. The diffraction limit for our imaging system was 
about 2 /im. For the determination of local quantities in 
the profiles, we averaged over ±5 jam around a given 
position. For temperature determination, the averag- 
ing region was restricted to an axial sector of ±60° to 
avoid corrections due to radial anharmonicities (see ap- 
pendix). The relative temperature is determined as 
V = T/Tfo = (-6L^3(-C))"^/^ where C is the fugac- 
ity obtained from the fit {Lis{z) = Y^^^i^^/k^ is the 
Poly-Logarithmic function of order s). 

Thermometry of ultracold Fermi gases 

In our previous work |T3 [19], temperatures have been 
determined by fitting the spatial wings of the majority 
component after expansion. However, we found that one 
can neglect collisions with the minority atoms in the core 
only for large population imbalances. In a simplified pic- 
ture, one can regard collisions with the inner core as col- 
lisions with a moving wall, which moves outward radi- 
ally and inward axially (due to the magnetic trapping 
potential). This results in different average kinetic ener- 
gies (transversely and axially) of the free majority atoms 
in the outer region. Figure [7| shows the density distri- 
bution of the majority and minority components after 
expansion. Although the temperature has been overes- 
timated by only 20% for typical experimental conditions 
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FIG. 7: Expansion of a population-imbalanced Fermi mix- 
ture. The absorption images of (a, d) the majority and (b, 
e) minority components were taken along (a, b) the axial z 
and (d, e) transverse y directions after expansion, (c) The 
azimuthally averaged column density profiles of the majority 
(red) and the minority (blue) cloud are obtained from (a) and 
(b), respectively. The excess majority atoms in the outer re- 
gion interact with the core during expansion. The contour 
lines of the outer part of the majority cloud (color inset) are 
not elliptical and have the shape of a horse-track. This shows 
that the minority cloud pushes the outer majority atoms in 
the transverse direction, which is also indicated by the hump 
of the majority profile at the edge of the minority cloud. The 
population imbalance was b ^ 55%. 



{S ^ 60%) in refs [T8l[T9], we do not regard this technique 
as well-calibrated absolute thermometry. 

One other concept for thermometry determines tem- 
perature as the derivative of entropy with energy. So 
far, this concept could be implemented only for bal- 
anced fermion mixtures with certain approximations, and 
due to the need of determining a derivative, could only 
be used to obtain temperatures averaged over a certain 
range [22] . 

Anharmonicity of the trapping potential 

For the determination of temperatures from the spatial in 
situ profiles it was necessary to address the anharmonic- 
ity of the trapping potential. Our trap is generated by a 
weakly focused (beam waist w ^ 125 jam) infrared Gaus- 
sian laser beam (wavelength 1064 nm) near the saddle 
point of a magnetic potential. The total trapping poten- 
tial is given as 



2p2 

V{p,z) = /7oexp( -) 



2 ^ 2^ 



(2) 



where 



- y'^ . Note that gravity has been compen- 



sated by a magnetic field gradient. The axial confinement 
comes mainly from the magnetic potential with oscilla- 
tion frequency of fz = 23 Hz. The transverse magnetic 
potential is anti-trapping and limits the trap depth as 



U 



1 



m{27rfpyw\l 




-)), 



(3) 



~\ 1 1 1 1 1 — 

30 40 50 60 70 80 
Averaging angle (9) 



FIG. 8: The temperature of the cloud was determined for 
various angles of the averaging sector. For a large angle, 
the large-x region is included in the averaged profile, resulting 
in a broadening of the spatial wings and consequently higher 
value of the fitted temperature. The red line shows the results 
of a simulation using the same parameters as the experiment 
(A = fp/fz — 6.15, Tfo = 1 /iK and the trap depth U/kB — 
2 /iK). 



where fp is the transverse oscillation frequency in the 
central harmonic region. When the trap depth is com- 
parable to the Fermi energy of a sample, the transverse 
anharmonicity will affect the shape of the cloud. Al- 
though in our experiments, the inner core and the outer 
cloud had the same aspect ratio as the trapping poten- 
tial, indicating the absence of anharmonic effects, anhar- 
monicities were not negligible in the spatial wings used 
to determine the temperature. 

This issue was addressed by adjusting the angular av- 
eraging region (Fig. |8|. Since the trapping potential is 
only anharmonic for large p, we could reduce the effect 
by decreasing the angle of the averaging sector around 
the z direction. Both the experimental data and an ex- 
act simulation for an ideal Fermi gas show that the fitted 
temperature remains almost constant up to a certain an- 
gle and then increases when the averaging sector includes 
more of the transversely outer region. In our tempera- 
ture determination, we chose the averaging sector to be 
±60° which was large enough to create low-noise profiles, 
but kept the effect of the anharmonicities to below 10%. 
The ID fit to angularly averaged profiles was computa- 
tionally more efficient than a 2D fit to a selected region of 
the image. In a 2D fit, one could also include anharmonic 
terms in the fitting function. 

Polarized superfluid at finite temperature 

When the two spin components have a chemical po- 
tential difference 2/i, the BCS-type superfluid has 
two branches of quasiparticles with excitation energies 
\/{ek -/i)2 + A2 ± h where Ck = h^P/2m. At finite 
temperature, the superfluid is polarized due to the large 
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thermal population of the lower branch compared to the 
upper branch. An interesting situation arises when h 
becomes larger than A, i.e. the lower branch has neg- 
ative energy quasiparticles, implying that even at zero 
temperature the superfluid state would have a finite po- 
larization. Our experiments show he < A ai very low 
temperature, suggesting that a polarized superfiuid state 
exists only at finite temperature. The breached-pair state 
with he > A ai zero temperature has been predicted in 
a stronger coupling region (on the BEG side of the Fesh- 
bach resonance) . Since A gradually decreases with higher 
temperature, it might be possible to have he > A at fi- 
nite temperature, at least in the weakly-interacting BCS 
limit where A smoothly approaches zero at a second or- 
der phase transition point. One interesting problem is 
identifying this gapless region of > A in the phase 
diagram for various coupling regimes. 
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